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Abstract 

The lowest ("vector") and next-lowest ("scalar") bound-state masses of the massive 
Schwinger model have been determined recently to a very high accuracy numerically 
on the lattice. Therefore, improved results for these bound-state masses from analytical 
calculations are of some interest. Here, we provide such improved results by employing 
both standard and renormal-ordered (fermion) mass perturbation theory, as well as a 
consistency condition between the two perturbative calculations. The resulting bound- 
state masses are in excellent agreement with the lattice results for small and intermediate 
fermion mass, and remain within 10% of the exact results even in the limit of very large 
fermion mass. 
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1 Introduction 

The massive Schwinger model is QED in 1+1 dimensions with one massive fermion species 
||1|]-|^. It is, one the one hand, the simplest non-trivial gauge field theory. On the other 
hand, it nevertheless shares some properties with more complicated gauge field theories 
like, e.g., QCD in 3+1 dimensions. The model has, for instance, a non-trivial vacuum 
|jl|, 0, and the fermions are confined in the sense that they do not show up as physical 
particles in the particle spectrum f^, |^ . 
The Lagrangian density of the model is 



Here, both the fermion mass m and the coupling constant g have the dimension of a 
mass, therefore the dimensionless parameter m/g effectively controls the behaviour of 
the theory. 

One interesting feature of the model is the formation of bound states. For small 
coupling (i.e., large m/g), the bound states are just fermion-antifermion bound states 
that are weakly bound by the 1+1 dimensional Couloumb potential [Q]. Therefore, the 
masses of the bound states are 2m plus some small (positive) binding energies and should 
approach 2m in the limit of vanishing coupling g. 

For large coupling, the model is best described in an equivalent bosonic formulation 
im, 1^. In fact, for a massless fermion (m = 0), the model may be solved exactly, and 
it is equivalent to the theory of one free, massive boson field ( "Schwinger boson" ) with 
Schwinger mass /io = 

g/V^ [0 i- For small but non-zero fermion mass (m ^ 0), this 
Schwinger boson turns into an interacting particle, and its mass acquires corrections 
which can be calculated perturbatively in m p, 0, |ll|. This Schwinger boson may be 
identified with the lowest bound state, the vector state. Further, a bound state of two 



Schwinger bosons is formed within the bosonic formulation of the theory 0, 11, 12|, and 
this bound state of two bosons may be identified with the next lowest bound state, the 
scalar state. 

These two lowest bound states remain stable for all values of the parameter m/g. 
In general, all higher bound states become unstable for sufficiently large coupling (small 
m/g) [0. For vanishing vacuum angle {6 = 0) there exists a third stable bound state 0, 
which may be interpreted as a bound state of one Schwinger boson and one two-boson 



bound state within a bosonic formulation |TT| . Here we will consider the case of vanishing 



vacuum angle (because all lattice data are for this case), but we will be concerned only 
with the two lowest bound states, because to our knowledge there are no lattice data 
available for the third bound state. 

For the calculation we will use perturbative expansions of the bosonic formulation of 
the theory for different choices of the normal-ordering mass, together with some consis- 
tency conditions among these different perturbation expansions. Qualitatively this has 



already been discussed in , therefore we will describe the calculation only briefly before 
presenting explicit results. 



2 



2 Vector mass calculation 



The theory of Eq. (|l]) is equivalent to the bosonic theory |l|, ||, || (for conventions we 
refer to [^) 



-N, 



A»0 



$(□ - H /io"^ (cos(V^$) - 1 

27r ^ 



(2) 



Here /iq = fi'/^r^/^ is the Schwinger mass of the massless (m = 0) Schwinger model, 
7 = 0.5772 is the Euler constant and N^j^ denotes normal-ordering w.r.t. /x. Further, a 
constant vacuum energy density has been subtracted, as usual. 
By making use of the well-known normal-ordering relation 



A'' 



(3) 



(see ) , we may rewrite the Lagrangian density (|^) for arbitrary normal-ordering mass 
/Z like 

-$(□ - H /im (^cos(v47r$) - 1 



-AT, 



(4) 



By shifting the quadratic (in $) part from the interaction part ~ cos A/47r$ to the free 
part of the Lagrangian (§), and by insisting that the normal-ordering mass is equal to 
the total mass of the free Lagrangian, we arrive at 



L = -Ni 



Mo 



-$(□ - H mfiQ cos(v47r$) + 27r$^ - 1 



27r 



(5) 



with 



Ho = e"'m + y fJ'O + e^'^w?. (6) 

The resulting Lagrangian (H) has a slightly more complicated interaction term (containing 
a $^ term in addition to the cos -\/47r$ term), and the dimensionless expansion parameter 
is now 

m 

(7) 



^ _ m 



This new expansion parameter A has the advantage that it tends to a finite and rather 
small value even for large mj g. 

The lowest bound-state mass M\ (the "vector" state mass or Schwinger mass) may be 
calculated from the two-point function ($(x)$(y)) for both Lagrangians (0) and (|^) by a 
perturbation expansion in the respective interaction terms. Via a Fourier transformation 
the mass pole may be determined. For ordinary mass perturbation theory (i.e., from the 
Lagrangian (^), the result up to second order is 



where 



A. 



(Mr')' = /i^ + 2e^ /iom + e^^m'^A^ + AS) , 

^^^L-ai^oM _ 1 ^ /o(r)(±e^'^o(") T 1 + 2iro(r)) 



(8) 
(9) 
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A+ = -0.6599, A_ = 1.7277 (10) 
see e.g. |P, 0, [ill]. The notation M™'^ means that the subscript 1 stands for the first 



(lowest) bound state mass, and the superscript m, 2 stands for the m perturbation ex- 
pansion up to second order. Here, /Xq is the order zero result, and there are both first 
and second order contributions. This result describes the vector state mass quite well for 
small fermion mass m, see Fig. 1, but it diverges, of course, in the limit of very large m. 

Starting from the Lagrangian (^, the order zero result is /xq. There is no first order 
contribution (because there is no $^ term in the interaction Lagrangian), and the total 
mass up to second order is 



(M^')2 =fil + e^^m\A+ + A_) - 2e^^ 



171^ 



/i^ + 2e^ cos^mA//i2 + e27 cos^ + e '''m (yl+ cos2^^ + y4_), (11) 



see 



13| . Please observe that one piece of the second order (in A) contribution, namely 



the term —2e'^'^m'^, has precisely cancelled the ©(m^) piece of the lowest order (in A) 
contribution, /Xq. This cancellation ensures that (Mj^'^)^, when Taylor expanded in m 
up to second order, coincides with the second order result of mass perturbation theory, 
(M™'^)2. (Mj^'^/m) reaches a finite value (which is about 55% above the true value 
2) in the limit of large m. Further, Mj^'^ describes the data quite well for small m/g, 
but already for m/g = 1 it is about 30% above the true value. Therefore, a further 
improvement would be desireable. This aim can be achieved by requiring that the result 
of the A perturbation expansion, when Taylor-expanded in m, should agree with the result 
from the m perturbation expansion also in higher orders. This requirement is perfectly 
sensible, because A has a well-defined Taylor expansion into m to all orders, and the 
physical results from both formulations of the theory should agree after all. A further 
Taylor expansion in m of the mass (Mj^'^)^ contains, e.g., the term 

fioliri^ = e , (12) 

in third order in m, and infinitely many higher order terms. Now the crucial point is 
that such terms cannot be produced by the perturbation theory in m, because Wick 
contractions of the interaction term cos v^47r$ all must contain integrals of a special 



type of transcendental functions (the E± functions of |TT|, |1^). Therefore, the higher 
orders in the A perturbation expansion must contain terms that precisely cancel terms 
like (p!^ ) when Taylor-expanded in m, as we have already observed for the second order 
contribution. The term (|12D is, e.g., precisely cancelled by the third order term (in A 
perturbation expansion) 

6{M^'')^ = -e^^'filX' = -e^^— (13) 

A* 

which must, therefore, be present in the A perturbation expansion. This cancellation has 
to continue in higher orders, and the Taylor coefficients of higher order terms, like 
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have to be cancelled, too. This cancellation condition enables us to perform a further 
partial higher order calculation in A perturbation expansion up to some order by 
employing the equation 

(M^^)2 ^ {Mty -filYl -nCiX' = {MT'^f + 0(m^+^) (14) 

where a Taylor expansion in m up to order is understood on the l.h.s. The coefficients 
Q can be determined recursively, which can be easily done with the help of Mathematica. 
We show the result for = 30 in Fig. 1 and in Table 1. As can be seen, the improvement 
for intermediate values of m/ g is dramatic. For m/g = 1, e.g., the result deviates from 
the lattice result by less than 0.5%. In the limit of very large m, the deviation from the 
exact result 2m is about 8%. 

Within A perturbation theory, the higher order terms that we have found from con- 
dition ( p^ correspond to disconnected products of Wick contractions, where in each 
connected component one cos$ term is contracted with one or several $^ terms of the 
interaction Lagrangian. There exist other types of terms which do not appear in the m 
perturbation expansion (contractions of two or more cos$ with at least one $^), but 
these may be detected in a way similar to the above argument only when compared 
with higher orders in m perturbation theory. Unfortunately, already a full third order 
calculation (both in m and in A perturbation theory) is rather involved. 



3 Scalar mass calculation 

The scalar mass may be calculated from the four-point function (nf=i ^{^i)) by focusing 
on the s-channel contribution and by again determining the mass pole. Here, we want 
to use the results of the last section for the vector mass, i.e., we want to use the masses 
(M™'^)^ and [M^'^Y for the respective perturbation theories, because they are better 
approximations to the physical masses than the bare masses fio resp. fio- However, naively 
using these masses within the Lagrangians (^ and (^) would lead to the problem that 
the propagating masses are no longer equal to the normal-ordering masses, which would 
produce tadpole contributions. Therefore, we change the normal-ordering masses in both 
Lagrangians and use for the subsequent calculations the Lagrangians 



-AT, 



Ml 



$(□ - M^)^ + —Mim (cos(A/47r$) 



(15) 



and 



-AT, 



Ml 



,2 ^ 



Ml)<^ + ^mMi (cos(y4^<l') + 27r<l>2 - l) + e^m(/xo - Mi)^"' 



^ . 6) 

where Mi will be finally chosen equal to Mf*' in the first case and equal to in the 
second case. Please observe that the $^ part of the interaction Lagrangian in ([l6| ) is not 
changed by the renormal-ordering and that we do not introduce a further shift such that 
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all <l>^ contributions belong to the free part of the Lagrangian. Here we just remark that 
this choice is perfectly legitimate, i.e., the resulting perturbation theory is as well-defined 
as the others. Later on it will become clear why this choice is even advisable for our 
purposes. 

Now, we have to calculate the bosonic four-point function, perform the Fourier trans- 
formation, focus on the s-channel contribution and find the mass pole in the kinematic 
variable s. The calculation is essentially equivalent to the calculations within mass per- 
turbation theory that have been done in [0, |ll|], therefore we just quote the result. One 
finally arrives at the pole mass equation (or gap equation) for the scalar mass M2 in 
leading order 

e'^m 4 1 1 

: arctan , (17) 



where y = M2 / Mi . It turns out that in leading order the gap equation is formally the same 
for both perturbation theories, the only difference being the different choices for Mi. The 
Ml masses have been determined fully up to second order in the respective perturbation 
expansions. It turns out that the above gap equation induces a leading order contribution 
to the mass ratio y which is of second order (see Eq. (|19]) below), therefore it is consistent 
to use the gap equation in leading order only. Further, we understand now why we had to 
choose the Lagrangian (|1^) as it stands. When Mi = M™'^ in (further renormal-ordered) 
m perturbation theory and Mi = Mi'^ in (further renormal-ordered) A perturbation 
theory, then it is automatically guaranteed that the resulting masses M2 will agree up 
to order when Taylor-expanded in m. Here we should be somewhat clearer about the 
Taylor expansion. The point is that we want to use the gap equation (|r^) directly for a 
numerical evaluation, which is then no longer accessible to a Taylor expansion. However, 
it is easy to generate an expansion in the coupling constant from the gap equation (|T^. 
Just rewrite it as 

,e''m 1 e 9,9 /.^x 

e = 4 arccot , = 4 - y'' (18) 

and solve this equation iteratively, starting with e = on the right-hand side. For the 
resulting expressions, we demand that they must agree up to order N when Taylor- 
expanded in m. 



There are several reasons for using the gap equation (17) directly for a numerical 
evaluation instead of a further expansion (iteration of Eq. (0)). Firstly, the leading 
second order contribution in m has a rather large coefficient ||12|, |lT], |15], 

f - 2 + 24.625^!^ (19) 

Ml 

and, therefore, approximates the lattice results rather badly already for quite small values 
of m. Secondly, after a further expansion, it is not even guaranteed that y < 2 (as it has 
to be for a bound state), whereas y < 2 holds automatically for the gap equation (|T7|) . 
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Remark: Please observe that for an iteration of Eq. (0) in leading order, as has been 
done in ^ , the further normal-ordering of Eq. ([TED is unnecessary, because it only 
changes results in higher orders. 

The gap equation has been evaluated numerically with the help of Mathematica. Al- 
ready within mass perturbation theory (i.e., for Mi = M™'^) the results of the numerical 
solution are quite reasonable up to m/g ~ 1/2 and are much better than the leading 
order result (JT^). For the Lagrangian (JT^), i.e., for Mi = Mj^'^, the results are shown in 
Table 1, again for the choice = 30. The results agree very well with the lattice data 
for small m. For m/g ~ 1, the discrepancy is about 4%, and in the limit of very large m 
the deviation from the exact result M2 = 2m is about 9%. 



4 Summary 

In this paper we have performed an improved calculation of the two lowest bound-state 
masses of the massive Schwinger model (the vector and the scalar bound-state mass). 
For this purpose we used the bosonized version of the theory with two different normal- 
ordering prescriptions and with two slightly different interaction terms in the Lagrangian, 
see Eqs. (0) and (||). In addition, we used a consistency condition between the two per- 
turbation theories in order to partially determine higher order contributins to the second 
(A) perturbation theory, see Eq. (|1^). The resulting masses agree very well with the re- 
sults from recent lattice calculations |jl5|, |l^ for small and intermediate values of the 
parameter m/g. In the limit of very large m/g they deviate from the known exact results 
by less than 10%. There are additional results on the masses Mi and M2 in the litera- 



ture, but we chose to compare with the lattice results of ||T^, [T^l as they seem to be the 
most accurate numerical estimates which are available up to now. For additional results 
on these masses with various methods, please consult [|15| and the literature cited there 
(some further references can also be found in [0). 

In principle, the results obtained in this paper could be further improved by methods 
that are analogous to those used in this paper. However, any further improvement requires 
a full third order calculation of the vector mass Mi within mass perturbation theory, 
which is quite difficult to do, as already mentioned. Further, the methods emlployed in this 
paper may, in principle, be used for the calculation of further observables in the massive 
Schwinger model, like, e.g., higher masses or decay widths and scattering cross sections. 
Also, it should be possible to use a variant of this renormal-ordered perturbation theory 
for a discussion of the multi-flavour Schwinger model where ordinary mass perturbation 
theory fails due to infrared divergencies |]T^ |1^ . 
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Ml 




m/ g 



Figure 1: The vector mass Mi, normalized such that the exact value is 1 in both limits 
m — i> 0, m — >• oo. The dots are the lattice results of |jl6|, the upper curve is the result M™'^ 
of mass perturbation theory, the lower curve is the improved result M^'^ with N = 30. 



Table 1: Comparison of vector and scalar masses Mi/g, M2/ g as functions of m/g. Our 

For m = 0, "this work" 



results are compared with the finite-lattice estimates of 
are the known exact results. 

m/g Mi/g Mi/g M^jg M^l g 

this work lattice, |T^ this work lattice, p!5| 





0.564 


0.564 


1.128 


1.11 


0.125 


0.789 


0.789 


1.470 


1.47 


0.25 


1.016 


1.019 


1.730 


1.74 


0.5 


1.474 


1.487 


2.165 


2.20 


1 


2.431 


2.444 


2.988 


3.12 


2 


4.486 


4.398 


4.876 


5.00 


4 


8.757 


8.340 


9.072 


8.85 


8 


17.40 


16.29 


17.78 


16.68 


16 


34.75 


32.24 


35.41 


32.56 


32 


69.46 


64.19 


70.71 


64.45 



